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We study solutions to the effective equations for the Bianchi IX class of spacetimes 
within loop quantum cosmology (LQC). We consider Bianchi IX models whose mat¬ 
ter content is a massless scalar field, by numerically solving the loop quantum cos¬ 
mology effective equations, with and without inverse triad corrections. The solutions 
are classified using certain geometrically motivated classical observables. We show 
that both effective theories -with lapse N=V and N=l- resolve the big bang sin¬ 
gularity and reproduce the classical dynamics far from the bounce. Moreover, due 
to the positive spatial curvature, there is an infinite number of bounces and recol¬ 
lapses. We study the limit of large field momentum and show that both effective 
theories reproduce the same dynamics, thus recovering general relativity. We imple¬ 
ment a procedure to identify amongst the Bianchi IX solutions, those that behave 
like k=0, 1 FLRW as well as Bianchi I, II, and VIIq models. The effective solutions 
exhibit Bianchi I phases with Bianchi II transitions and also Bianchi VIIq phases, 
which had not been studied before. We comment on the possible implications of 
these results for a quantum modification to the classical BKL behaviour. 

PACS numbers: 04.60.Pp, 04.60.Bc, 98.80.Qc 


I. INTRODUCTION 

In recent years various homogeneous cosmological models have been studied within the 
context of loop quantum cosmology (LQC) [1]. In particular, for the models that have 
been exactly solved at the quantum level, the dynamics of sharply peaked, semiclassical 
states is very well described by an effective theory that incorporates the main quantum 
corrections to the dynamics [2H5]. These examples give one confidence that the effective 
LQC description will provide reliable information about semiclassical states for models which 
have not been fully solved. It is in this context that we study the Bianchi IX model. Within 
homogeneous and anisotropic models of the universe, the Bianchi IX model exhibits very 
interesting features. One of them is due to the positive curvature of the underlying spatial 
manifold, which has the consequence that all classical solutions exhibit a recollapse for 
generic matter content following Einstein’s dynamics. At the effective level that we shall 
consider here, an inhnite number of bounces and recollapses appear, a feature that was 
hrst seen in the Friedman-Lemaitre-Robertson-Walker (FLRW) model with k = 1 [HHTn]. 
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Furthermore, in the Bianchi IX model, the so called inverse triad corrections of LQC play 
an important role and modify the quantum theory yielding different descriptions. A valid 
question is whether there exist physical criteria to select one theory over another, as was done 
in the isotropic case In other words, can the semiclassical limit choose the ‘physically 

correct’ theory? After all, the semiclassical limit of the quantum theory is expected to make 
contact with the classical theory, and this can be useful to select the theory that satishes 
some physically motivated requirements. If we assume that the semiclassical limit of each 
quantum theory is well described by their effective theory, we can try to answer this question 
by studying the effective theory. One should note that this issue has already been explored 
for some models m |15]. Here we want to explore this issue further, and consider an 
improved quantization of the Bianchi IX model, put forward in [16], that incorporates the 
so-called inverse triad corrections. 

The answer that we hnd in the numerical exploration of the different effective theories, is 
that each of them describes different dynamics but, as we show, in the limit of large held 
momentum the dynamics of all the different theories is almost the same. For the Bianchi IX 
models, large held momentum amounts to spacetimes that attain a large volume (in Planck 
units) before re-collapsing. From this point of view, even if there are diherent quantum 
theories with diherent semiclassical limits, if the universe grows to be ‘large’ then both 
theories describe the same physics in the large held momentum regime. 

The Bianchi IX model has long been studied in the context of Loop Quantum Cosmology. 
The hrst study [T7|, was within the so called po dynamics, whose infrared limit has been 
shown to be unphysical mm- In order to solve the problems with the infrared limit, the 
improved fl dynamics for Bianchi IX was constructed ng. It is within the framework of the 
improved /2 dynamics that we consider the ehective theories for this work. We will show how 
to recover the correct infrared limit and resolve the big bang singularity in the ultraviolet 
limit. 

Another important feature of the Bianchi IX model comes from the BKL (Belinskii- 
Khalatnikov-Lifshitz) conjecture [TOH^ . which suggest that the local dynamics of inho¬ 
mogeneous cosmologies near the big bang singularity is dominated by the time derivatives, 
which become more important than the spatial derivatives. This fact makes the local dynam¬ 
ics at each spatial point independent and, therefore, can be approximated by a homogeneous 
universe. In general, the dynamics is described by a Bianchi IX universe, which exhibits an 
oscillatory behaviour between different Bianchi I phases with Bianchi II transitions. Another 
important feature of the BKL conjecture and the evolution near the big bang singularity is 
that, in general, the dynamical contribution from matter is negligible and the universe can 
be described as a vacuum universe. If one puts together these two elements that follow from 
the BKL conjecture, then one can conclude that the dynamics of any universe near the ini¬ 
tial singularity can be described locally as a homogeneous universe where all the dynamical 
contribution comes from the anisotropies or, in other words, a vacuum Bianchi IX universe. 
An important caveat is that this is true only if the matter content is different from a massless 
scalar held. In this particular case the dynamical contribution from the matter turns out 
to be important |22l |23|. This peculiar and unique behaviour of the massless scalar held 
justihes its study within the LQC Bianchi IX model. First, if the model under consideration 
consists of a scalar held together with another kind of matter, it is only the mass-less scalar 
held that contributes to the dynamics in a signihcant way. Second, if there is no scalar held 
in the universe, then the conjecture states that the dynamics is dominated by a vacuum 
universe, which can be seen itself as a limiting case of a universe with a massless scalar held. 
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when the momentum of the held goes to zero (for more details on the classical behaviour of 
Bianchi IX see |23l [2l])- As this discussion suggests, the massless scalar held together with 
the vacuum case, represent the most interesting cases to be studied in a Bianchi IX model. 
In this manuscript we restrict ourselves to the massless scalar case within LQC. We shall 
only briehy comment on the vacuum case. 


Since the Bianchi I and II models play an important role in the classical behaviour of Bianchi 
IX, we shall also consider the ehective models of Bianchi type I and II, highlighting the most 
important aspects of these models, with the purpose of understanding the Bianchi IX effec¬ 
tive dynamics. Additionally, given that the quantization process can produce unexpected 
results, we include also the study of the isotropic FLRW models with k = 0,1, in order 
to illustrate some new results in the dynamics of Bianchi IX. Finally, since our study is 
numeric, we use some known results of the isotropic models and Bianchi I |2Sl I2S], in order 
to test the validity of the numerical implementation. It is important to point out that these 
cosmological models are not sub-cases within the quantum Bianchi IX model [TS], but are 
limiting cases instead. A special case is given by the classical isotropic k = 1 model which 
is embedded in the Bianchi IX universe. For the LQC effective models, this is also the 
case only when the inverse triad corrections are not included. If one includes them, then 
the isotropic model is only a limiting case [ZIEIE]. Here we shall show that in, the large 
momentum regime, the isotropic limits of the effective Bianchi IX and the effective FLRW 
k = 1 coincide, regardless of whether there are inverse triad corrections or not. 


This paper is the second in a series. In the first one [16], different version of the loop 
quantization of Bianchi IX model are introduced. In contrast to previous treatments, the 
inverse corrections are fully implemented. The effective theory is constructed and some of 
its properties are analysed. In the third manuscript of the series [2Z| we shall present some 
discussion regarding qualitative features of the effective dynamics for generic Bianchi IX 
LQC models, including the vacuum case. One should also note that a summary of some of 
these results has already been reported in [28] . 


This article is organized as follows: In Sec. |TT|we review the relevant results for the classical 
Bianchi IX and discuss some classical observables relevant for the effective theory. Next, in 
Sec. [ITT] we summarize the two effective theories used in our study. In Sec. |IV]we discuss the 
numerical results of the effective theories. Since we do not have an analytical understanding 
of the space of solutions, nor of the “generic” behaviour, we explore several limiting cases. 
Finally, the conclusions are outlined in Sec.jVj Additionally, we include an Appendix with all 
the Hamiltonians and equations of motion used in this manuscript, and a second Appendix 
where we display the numerical convergence tests of our integrations. 


II. CLASSICAL DYNAMICS 

In this section we shall recall the classical description of the Bianchi IX model. This section 
has two parts. In the hrst one we review the basic variables used in both the geometric and 
the connection dynamics formulations. In the second part, we introduce the geometrical 
invariants and observables that are useful to extract physical information. 
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A. The model 


In the Bianchi IX model (like in the FLRW with k = 1) the spatial manifold S has the 
topology of a three sphere The physical metric on S is qab ■ = where the physical 

forms are and the hducial ones are [IB] 

°u:l = sin/3sin7(da)a + cos7(d/?)a , 

°ojI = - sin/?cos7(da)a + sin7(d/5)a , 

°u:l = cos/3(da)a + (d 7 )a , 


with a G [0,27r), [3 G [0, tt) and 7 G [ 0 , 47 r). The radius of the three sphere as measured 
by the hducial metric is oq, therefore the volume with respect to the hducial metric is 
Vo = 271^ al- We dehne io := , that can be rewritten as lo ='■ 'dao, with := (27r^)^/^. 

The triads and connections are 


4* = — 1® 
to 




^0 


o a 

■ 


Using these dehnitions and the lapse function N = V = ^ypiP 2 P 3 (with pj > 0), the Hamil¬ 
tonian constraint is [ 18 ] 


nCla. 


1 


P1C1P2C2 + P2C2P3C3 + P3C3P1C1 + 2d[piP2C3 -f P2P3C1 + P3P1C2) 


+^9^(1 + 7 ') 

- 0 . 



Poisson brackets of the basic variables are {c®,pj} = SttGjSj and = 1. The equa¬ 

tions of motion come from the Poisson bracket of the basic variables with the Hamiltonian 
constraint 


Pi {Pi^ ^Hbix} 5 ^ ! 

0 ^Hbix } ’ ^Hbix } ’ 


which are given by 


P<p = 

0 = 

Pi 


Pi = 


7 


P2C2 + P3C3 -1- 2'd 


P 2 P 3 

Pi 


Cl = — 

7 


P2C1C2 + P3C1C3 2'd(p2C3 -F P3C2) 

V m m 2p| j 


The other equations can be obtained by permutations of the labels. 
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B. Observables 


In order to study the homogeneous cosmological models, there are two interesting physical 
time variables. The hrst one is harmonic time r, which corresponds to the choice of lapse 
equal to the volume, N = V. The other one is the cosmic (proper) time t (corresponding to 
the choice of lapse iV = 1). These two time variables are related by the equation 


d _ 1 d 

-Jp-iPm dr' 


( 2 . 1 ) 


The observables are usually dehned with respect to the cosmic time. The derivative with 
respect to the cosmic time is denoted as O' = dO/df, therefore O' = 0 /a/PiP 2 P 3 , where the 
‘dot’ stands for derivative with respect to harmonic time r (we are using the notation from 
references [HI [32l [33] )• In order to understand the singularity resolution and the evolution 
of the system as dehned by the classical and effective equations, we study the following 
observables: 


1. Directional scale factors, 




PjPk 


0 V Pi 


which come from the relation pi = ajOk^Q, with i ^ j ^ k ^ i and Pi, Oi > 0. 
2. Directional Hubble parameters, 


with i ^ j ^ k ^ i. 
3. Expansion 


Oi 2 \pj pk Pi 


e = y = H, + H, + H3. 


An equivalent function is the Average Hubble parameter 

H=^- = ^{H, + H2 + Hs). 


4. Matter density 

pI - pI 

21/2 2piP2P3 ' 

The dynamical contribution due to the matter is measured by the Density parameter 


SttG p 

~^1p' 


( 2 . 2 ) 


5. Shear 


= -[(TTi - H2Y + - H3) 


+ {H-i - H^y-] A 


2 = 1 


^ Note that this dehnition of H differs from the standard definition H = ^(Tahcr®^. 
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The dynamical contribution due to the anisotropies is measured by the Shear parameter 


: = 


3(7^ 




(2J) 


6. Rieei sealar 



( 2 . 4 ) 


The p' and 6 ' are computed from the equations of motion. In the classical case 9' is 
given by the Raychaudhuri equation 9' = —\9‘^ — — WnGp. 


7. Intrinsie eurvature, one feature of Bianchi II and IX models is that the spatial curva¬ 
ture is different from zero. The intrinsic spatial curvature is given by 


where 


P2P3 

Xi — O'IA ^ , X2 — 0X2 

pi 


- 2(a;ia;2 

X1X3 + X2X3)] , 

( 2 . 5 ) 

IpiPs 

P 1 P 2 

(2.6) 

V pi 

^3 — Oi3\ 3 • 

V P 3 


The values of a,- are: 


• Bianchi I, ai = 0(2 = «3 = 0. 

• Bianchi II, ai = 1, 0:2 = 0, 0:3 = 0 or permutations. 

• Bianchi IX, cki = 02 = 0:3 = /q, with /q = 

• The isotropic FLRW model with k = 1, is obtained from the Bianchi IX model 
by setting xi = X 2 = X 3 . 


One can also dehne another quantity that gives information about the dynamical 
contribution of the intrinsic curvature, named Curvature parameter, 


K 


202 



2 I 2 

X2 + X3 


2{xiX2 + X1X3 + a;2a;3)] , 


( 2 . 7 ) 


The parameters O, T? and K satished the classical relation 

VL + T? + K = 1 . ( 2 . 8 ) 

As we shall see later, these parameters are ill dehned at the quantum bounce (where 
0 = 0). Therefore they are only useful away from the bounce, and in the classical 
region. 

8. Kasner exponents 

k. = ( 2 . 9 ) 

These parameters are useful to determine when the solutions are of type Bianchi I. The 
sign of Hi determines whether the direction is expanding [ki > 0) or contracting 
{ki < 0 ). 
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The classical Bianchi IX model does not have complete exact solutions, so there are only 
some general results. We summarize those results that are important to us (a more complete 
review can be found in [23]): 

• All the Bianchi IX solutions recollapse if the matter satishes the dominant energy 
condition and has a non-negative average pressure [ 22 ] • 

• Bianchi IX has, as a limiting case, the Bianchi I model when Xi —>■ 0, with i = 1, 2, 3. 

• Bianchi IX has, as a limiting case, the Bianchi II model when xi 7 ^ 0,^2 —t 0 , 0:3 —?• 0, 

and permutations of xi,X 2 ,x^. 

• Bianchi IX reduces to fc = 1 FLRW when pi = P 2 = Ps and Ci = C 2 = C 3 . 

• Bianchi IX has, as a limiting case, the k = 0 FLRW when pi = P 2 = Ps, ci = C 2 = C 3 

and Xi —)■ 0, with i = 1, 2, 3. 

• The solutions to the vacuum Bianchi IX can be approximated near the big bang 
singularity as evolving in Bianchi I phases with Bianchi II transitions |24j . 

• Bianchi IX reduces to Bianchi VIIq when ai = 0, 0:2 = 1, 0:3 = 1 (and permutations) 

and is a limiting case when ai = 0:2 = 0:3 = one of the Xi goes to zero, with 

i = 1,2,3. 

The Xi variables will be very useful when we study the evolution of the effective solutions of 
Bianchi IX, because they give information about the transitions of the solutions. In order 
to distinguish between the transitions that Bianchi IX has, we chose Bianchi IX {ai = a 2 = 
as = /q) and study the evolution of the Xi variables. The different limits that the Bianchi 
IX solutions can approach may be classihed as follows, 

• Bianchi I: {xi —)■ 0 ,a :2 —)■ 0 ,X 3 —?• 0). 

• Bianchi II; (xi —0, 0:2 —)• 0) or (xi —)■ 0, 2:3 —j- 0) or {x 2 —)■ 0, X 3 —0). 

• Bianchi VIIq: (xi —)■ 0 ) or (x 2 —)■ 0 ) or (X 3 —)■ 0 ). 

Note that Bianchi VIIq is included in this classihcation given that, classically, it forms part 
of the Bianchi IX dynamics. At the effective level it is not known whether this is still 
true, therefore it is good to include this case in our study. It is important to emphasize 
that, so far, there does not exist a Loop Quantum Cosmology version of Bianchi VIIq and 
therefore neither its effective version. Then, when we describe a particular regime behaving 
as Bianchi VIIq, this refers to solutions in which one of the x, goes to zero. It would 
probably be worthwhile to quantize the Bianchi VIIq model, obtain its effective dynamics, 
and compare with what we get here for Bianchi IX and call Bianchi VIIq. 

The evolution of each Xj can give us information on how close to other Bianchi models is 
the solution of the full Bianchi IX dynamics. Moreover, the values of hi, and K can give 
us information about how the dynamics of the system is affected due to the contributions 
from matter, anisotropies and curvature. Classically, the identity hi -|- -|- iC = 1 is always 

satished, and this can be used to verify when the evolution of the system is approaching the 
classical region. Furthermore, the Kasner exponents fy can be helpful to classify which kind 
of Bianchi I solutions we get. 
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All of these elements will be useful to study the solutions to the effective equations. Even 
when our goal is to study the Bianchi IX dynamics, due to its complexity it is necessary to 
study the simpler Bianchi models as well as the isotropic models. The effective Hamiltonian 
of each of these theories are shown in the next section. All the equations of motion are 
displayed in Appendix 


III. EFFECTIVE DYNAMICS 


The “improved dynamics” quantization of the Bianchi IX model [IH] consists in calculating 
the connection from the holonomies and use this connection to dehne the curvature, which 
is then promoted to a quantum operator. There exists a freedom in the choice of the lapse 
function; this freedom is not important physically at classical level, but at quantum level 
it dehnes different quantum theories. In this manuscript we consider two effective theories 
that come from two different quantum theories for the Bianchi IX model. The hrst choice 
is when the lapse function is equal to the volume N = V, and the second one when N = 1. 
The effective Hamiltonian for N = V is IT8l 


^(1) __ 


P 1 P 2 PZ 


sin jjLiCi sin p, 2 C 2 + sin ^202 sin P 3 C 3 + sin ^303 sin piC\) 


-' -sm/i3C3H- ^z^sm^iCi-\ - ^^^smfi2C2 I 


V y/Ps 
^9^(1 + 72 ) 


7-., 


2{pi+pt+pi) 


y/P~l 

P1P2 

P3 


y/n 

P2P3 \ ^ f PsPl ^ ^ 




Pi 


P2 


(3.1) 


with'd = ( 27 r^)^/^, = d-y/STry^pj, 7 the Barbero-Immirzi parameter and 

fii = ——) p2 = — —) fis = 


P3 


P2P3 V P 1 P 3 V P 1 P 2 

When the lapse function is iV = 1 the effective Hamiltonian is given by dniEO] 


C 


( 2 ) _ V^A{V)h\V) 


sin piCi sin ^202 + sin piCi sin ^303 + sin ^202 sin 


87rGI/672A2 

dA { y ) h \ V ) / 2 2 . 22. 22- ^ 

■ /^^QyA^2x I ^ 1^2 Sin P3C3 + P2P3 sin PiA + P 1 P 3 sm P2C2 J 

i9^(l + j^)A(V)h^(V) 


SttGV/X'^ 
2VIp^, + pI + pI] 


X 


(piP2y + (piPs^ + {P 2 P 3 Y 


h%V)\ 


h\y)y\ 2^. 

P 4 , ~ 


21/« 


0 . 


(3.2) 
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with Vc = 27r7A£|[ and 


h(V) = Vv + v,-^\v-v,\, 

A(V) = ^jV + V,-\V-V,\). 

The Poisson brackets are {c\pj} = SnG'jSj and {0,p</,} = 1. The equations of motion are 
shown in appendix [A| 

Given that the classical Bianchi IX model has phases of Bianchi I and transitions of Bianchi 
II [21], it is interesting to study what happens at the effective level. Furthermore, we also 
consider the effective Hamiltonians for Bianchi I |3T] and Bianchi II [32], which can be 
written as [32] 


Chbii - 


P 1 P 2 P 3 

87rG72A2 

1 


sin /iiCi sin /i 2 C 2 + sin P 2 C 2 sin P 3 C 3 + sin P 3 C 3 sin piCi 


SirG'-y'^ 


a{p2P3f^^ . _ n , 2 ^ ap2P3\ 

— sm/iiCi - (1 + 7 ) —- 




V 2pi J 


-d«o, 

2 


(3.3) 


where a plays the role of a switch between Bianchi I (a = 0) and Bianchi II (a = 1). The 
solutions to the effective equations derived from this Hamiltonian were already studied in 
[3H] . and we shall use them to study the solutions of the effective Bianchi IX models. 


The two effective Hamiltonians of Bianchi IX ( 3.1[3.2 ) can be seen as the generalization of 
the two effective theories for the closed FLRW {k = 1), that are obtained using the same 
quantization method for the curvature and the two lapses N = V and X = 1 [3 IE]- The 
effective Hamiltonian for the closed FLRW model with N = V is given by [ 6 ], 


C 


( 1 ) 

Hk=i 


31/2 


[ sin2 A/3 - 2T) sin A/3 + (1 + 



(3.4) 


with D = X'd/V The effective Hamiltonian with iV = 1 is given by [H], 


C 


( 2 ) 

Hk=i 


3A{V)V 

SttGj^X'^ 


[ sin^ A/3 — 2D sin A/3 + (1 + 



(3.5) 


The quantum theory from which this effective Hamiltonian comes includes more negative 
powers of the volume, i.e., inverse triad corrections. A more detailed discussion about these 
quantizations and the solutions to their effective equations can be found in [T] El EO] . 

The variables (V, /3) are related with (c, p) by [3]: 


V = p^/\ /3 = ^. (3.6) 

Vp 

The Poisson brackets are 

{/3,1/} = dvrGy , {(/),p^} = l. (3.7) 

Finally, we want to consider the effective Hamiltonian for the flat FLRW model k = 0, 
since it is a limit of Bianchi IX spacetime and can be useful to check the correct numerical 
implementation. Furthermore, it can also be used to understand the behaviour of Bianchi 
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IX near to the isotropic universes. The effective Hamiltonian for FLRW with k 
by 






0 is given 


(3.8) 


with lapse N = V, and V the physical volume of the fiducial cell. 

Once one gets all the effective theories, the next step is to study their numerical solutions 
and use them to classify and understand the Bianchi IX dynamics. This will be reported in 
the next section. 


IV. NUMERICAL SOLUTIONS 


In this section we study the numerical solutions to the effective equations of the LQC Bianchi 
IX model. We consider two effective theories of Bianchi IX, one with lapse N = V [IH] , and 
another with iV = 1 (TB]. Given the difficulty to study the solutions of these theories our 
strategy follows two paths. The first one is the discussion of the two effective theories from 
the analytical point of view. Second, we study the space of solutions where the two theories 
describe the same physics and justify the relevance of this region in solution space. Next, we 
consider several regimes that the solutions space may have, in order to exhibit different types 
of behaviour. We do not expect those cases to be exhaustive, but they give us important, 
physical, information. In particular, we consider the Bianchi I limit, the non-shear limit and 
the isotropic limit, within this set of solutions. 


One of the objectives of this part is to answer some questions regarding the effective solutions. 
First, we would like to know whether the singularity is resolved and whether the solutions 
reduce to the classical ones far from the bounce. As we shall describe in detail, the answer 
is in the affirmative. Second, we would like to verify that Bianchi IX reduces to the closed 
FLRW {k = 1) when pi = P 2 = Ps and Ci = C 2 = C 3 . Third, we study the limit of flat FLRW 
{k = 0) when pi = P 2 = Ps, Ci = C 2 = C 3 and Xi, X 2 , X 3 —)■ 0. Finally, we investigate whether 
the LQC Bianchi IX model has (LQC) Bianchi I as a limit when xi,X 2 ,X 3 —?• 0. 

In order to answer the second question we must discuss the differences between the two 
effective theories from the analytical point of view. It can be shown [30] that the Hamiltonian 
Cy ^ for Bianchi IX with lapse N = V, Eq. (3.1), reduces to the effective closed FLRW with 

( 1 ) _ ^( 1 ) 


N = V, Eq. (3.4), when pi = P 2 = P 3 and Ci = C 2 = C 3 , in which case Chbix ~ 


Hk 


This 

means that, when one quantizes both the isotropic and the anisotropic models, the solutions 
to the anisotropic effective model must contain the solutions of the isotropic effective model. 
For the full quantum theory, one expects that the quantization of a less symmetrical model 
reduces to the quantization a more symmetrical one, as is explicitly the case for the Bianchi 
I model [26]. In the case of Bianchi IX such reduction is still to be shown at the quantum 
level. On the other hand, the new quantization for Bianchi IX with N=1 and inverse 
corrections put forward in [3 [30], has a different behaviour. When we consider its effective 


(2) A -k 

Hamiltonian Chbix’ ^0- (3.2), and choose pi = P2 = Pa, ci = C2 = C3, it does not reduce to 

the closed FLRW with N = 1, Eq. (3.5), namely Chbix ^ Therefore, the dynamics 

that describe the two Hamiltonians for the isotropic case are different. In this case the 
process of quantization and the symmetry reduction do not commute. A more detailed 
study regarding the differences of these two theories at the quantum and effective levels can 
be found in 
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As the detailed study of the closed FLRW model with different quantizations miH] has 
shown, their effective theories, while providing different dynamics for the bounce, yield the 
same results in the limit of large volume or large momentum of the held. The natural 
question is whether such feature extends to the anisotropic scenario. With this question 
in mind is that we study the numerical solutions of Bianchi IX. The idea is to look for 
solutions in the limit of large momentum of the held, i.e., where the volume at the bounce 
is expected to be much larger than the Planck volume (as is the case for the /c=0,l models). 
It is important to recall that in both the closed FLRW and the Bianchi IX model it is valid 
to talk about large volumes, since the spatial manifold is compact and has a hnite physical 
volume. This does not happen in neither of the Bianchi I, II and open FLRW models where 
the quantity V represents the physical volume of the hducial cell. 

This section contains hve parts. In the hrst one, we study the limit of large held momentum 
and show that the two effective theories of Bianchi IX reproduce the same dynamics. In the 
second part, we show that the effective solutions resolve the big bang singularity, and repro¬ 
duce the classical solutions far from the bounce. In the third one, we explore the isotropic 
limits, for both hat and closed FLRW models. This allows us to compare the numerical 
solutions with the known analytical and numerical results for the isotropic models. In the 
fourth part, we consider the limit of small shear, which shows a very diherent behaviour than 
the isotropic limit. In the hnal part, we study the Bianchi I limit, explore the transitions and 
discuss new results concerning to Bianchi IX and Bianchi VIIq. Throughout the units used 
to perform the simulations are: c = 1, h = 1, G = 1 ,7 = 0.23753295796592. The numerical 
method used to integrate the equations is a Runge-Kutta 4 method.^ All the plots show 
simulations done with the cosmic time (lapse N = 1). When we choose t = 0 as the time 
that the bounce occurs, we should clarify that this is the bounce in the geometric mean of 
the scale factors, or equivalently when the expansion 6 * (or average Hubble rate) vanishes. 


A. Large Field Momentum Limit 

First we are going to study the region in the space of solutions where both effective equations 
give the same dynamics. This region correspond to the solutions with large held momentum. 
In this limit we can say that the two theories are equivalent, in the sense that, given an 
interval of time tf — ti < 00 and a tolerance e > 0 , it is always possible to hnd initial 
conditions such that the solutions of both theories evolve within the tolerance e > 0 for the 
interval of time tf — ti < 00 . It is important to clarify that all the assertions about the 
effective solutions apply to the numerical solutions that we study, and not for all the space 
of solutions. Nevertheless, we expect that the qualitative results found in our numerical 
explorations can be extended to other, generic within the respective class, initial data. 

In Figs. mil the comparison between the solutions of both effective theories can be 
appreciated: the three scale factors have the same behavior in both theories. The bigger the 
volume at bounce is (or the momentum of the held), the closer the two effective solutions 
are. The initial conditions are: fiiCi = 7 r/ 3 , fi 2 C 2 = vr/2, paCs = 7 r/ 4 , pi = 13000, p 2 = 27000 
and Pa = 42000. The momentum of the held is p^ = 3.17 x 10®, that comes from the 
Hamiltonian constraint. 


^ The program is available up on request. 
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FIG. 1: Time evolution of the volume for the two different effective theories, with lapses N=V and 
A^=l. The initial conditions are: fiici = tt/S, fi 2 C 2 = vr/2, / 23 C 3 = 7 r/ 4 , pi = 13000, p 2 = 27000 
and p 3 = 42000. As it can be seen from the left hgure, the universe undergoes a series of bounces 
and recollapses, and the dynamics of both theories is indistinguishable. In the right figure we zoom 
into one of the bounces where the dynamics of both theories coincides again. 


Scale Faclors 



time 


Scale Factors 



time 


FIG. 2: Time evolution of the directional scale factors ( 01 , 02 , 03 ) for the two different effective 
theories, with lapses N=V and A^=l. The initial conditions are: flici = 7 r/ 3 , fl 2 C 2 = 7 r/ 2 , 
AI 3 C 3 = 7 r/ 4 , Pi = 13000, p 2 = 27000 and ps = 42000. Note that the dynamics of both theories 
coincides, even through the bounce (right). 


From now on, the study of the effective solutions will be at the regime of large held momen¬ 
tum. Since we have shown that in that limit both effective theories are indistinguishable, all 
the conclusions that we reach will apply for both effective theories of the Bianchi IX model. 
In particular, the big bang is replaced by a bounce in all cases. This will be the subject of 
study in the next section, together with the comparison with the classical solutions. 
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Density Shear 




FIG. 3; Time evolution of the density and shear for the two different effective theories, with lapses 
N=V and A^=l. The initial conditions are: fiici = vr/S, fi 2 C 2 = fisCs = vr/d, pi = 13000, 
P 2 = 27000 and p^ = 42000. Again, both theories can not be distinguished, even in the Planck 
regime. 




FIG. 4: Time evolution of the volume. Here we compare effective solution with the classical 
solutions they approach, before and after the bounce (that occurs at t = 0). It can be appreciated 
that the classical and effective solutions have the same behaviour in the classical region. Near 
to the bounce the classical solutions go to zero volume (big bang singularity) and the effective 
solution bounces. The initial conditions are; fiici = 37r/8, p, 2 C 2 = '7r/2, = Svr/S, pi = 15000, 

P 2 = 10000, p 3 = 20000. The initial conditions for the classical equations are taken from the 
effective evolution at the maximal volume, and evolved backwards. 


B. Classical Limit 

We want to verify that the classical big bang singularity is resolved by the effective dynamics 
and check that the classical dynamics is recovered in some interval of time. This add self- 
consistency to the effective theory, in the sense that it reproduces the classical dynamics far 
from the bounce (consistent with general relativity), and it solves the big bang singularity 
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Scale Factors 



FIG. 5; Time evolution of the scale factors (01,02,03). As in the previous hgure, we compare the 
effective and classical solutions. It can be appreciated that each direction bounces (at a different 
time) in the effective solution, while all the classical solutions go to zero at some point, either before 
and after the bounce (that occurs at t = 0). The initial conditions are: JliCi = Stt/S, p, 2 C 2 = vr/2, 
AI 3 C 3 = Stt/S, Pi = 15000, p 2 = 10000, p^ = 20000. The initial conditions for the classical equations 
are taken from the effective evolution at the maximal volume. 




time 


time 


FIG. 6: Time evolution of the density and shear. The effective and classical solutions are compared. 
It can be appreciated that, for the classical solutions the density and shear diverge, either as the 
solution approaches a big crunch singularity if it is contracting, or possesses a big bang singularity if 
it is expanding. Both scalar quantities remain finite in the effective solution. The initial conditions 
for the effective equations are: pici = Svr/S, p, 2 C 2 = 7r/2, = 57r/8, pi = 15000, p 2 = 10000, 

P 3 = 20000. The initial conditions for the classical equations are taken from the effective evolution 
at the maximal volume. 
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incorporating then what we expect to be the dynamics of the semiclassical states of the 
quantum theory. 

We choose the initial conditions near to the bounce for the effective equations and evolve 
them back and forward in time. The initial conditions are: JiiCi = Svr/S, J12C2 = 7r/2, 
PaCs = Stt/S, Pi = 15000, p 2 = 10000, pa = 20000. In order to do an easy comparison with 
the classical solutions, we take the values of the variables (cj,pj) at the maximal volume, 
that come from the effective evolution, and later on we introduce them as initial data for 
the classical equations, and then we evolve back and forward in time. 

The effective and classical solutions are compared in Fig. Here we show the time 

evolution of the total volume and the scale factors. Note that the effective and classical 
solutions coincide up to a point close to the bounce, then the classical solutions have a 
singularity, while the effective solutions bounce and connect with the classical solutions 
before and after the bounce. The recollapse is due to the positive spatial curvature and that 
holds for any matter content that satishes the dominant energy condition. The bounce is 
due to the loop quantum cosmology effects. Therefore the universe continue bouncing and 
recollapsing forever. 

Once we verify that the effective solutions reproduce the classical solutions and solve the 
big bang singularity, the next step is to study the space of effective solutions. Since the 
dynamics of Bianchi IX is non trivial and there are no analytical solutions for the classical 
dynamics, one is not expected to get a complete and fully detailed study of the effective 
solutions. Still, it is possible to perform a qualitative study and get a better control on the 
space of solutions. For this reason we study different limits of the solutions. One interesting 
limit is when the solutions are isotropic, which is the case we shall explore next. 

C. Isotropic Limit 

Within the Bianchi IX solutions, the ones that are close to the isotropic solutions are very 
important, since it is expected that any ‘realistic’ model of the universe must contain these 
solutions. After all, our universe is almost homogeneous and isotropic at large scales. Since 
the numerical solutions for the full quantum isotropic universes are known [21 El El El E5], 
they can be compared against the effective Bianchi IX solutions in the isotropic limit to test 
its validity. 

In hgurej^we show the evolution of the Kasner exponents ki,k 2 , k^ and their sum ki + k 2 + k 3 . 
It can be seen from the hgure that all the exponents have the same value (fc* = 1/3) and 
the total sum is one. This proves that the expansion rate in each direction is the same. 
The sign of the Kasner exponents indicates that its corresponding direction is expanding 
{ki > 0) or contracting {ki < 0), while the sign of the sum ±1 tells us the sign of the total 
expansion 6. The bounce or the recollapse is when 0 = 0, if 0 > 0 the universe is expanding 
and for 6* < 0 it is contracting. This dynamics can be seen in hgure [^ where it is shown the 
evolution of the volume and the expansion 6. Note that there are two types of bounces. The 
volume at the bounce differs for each type of bounce and it has two different values. This 
feature was hrst studied, for the k=l FLRW model in [7], and correspond to the dynamics 
of the effective theory associated with a connection based quantization (as opposed to the 
curvature based quantization of m)- This allows us to verify the results for the isotropic 
k = 1 case, and be conhdent of our numerical results. The initial conditions are: fiiCi = 7 t/2 , 
Pi = 100, with z = 1, 2,3. 
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time 

FIG. 7: Evolution of the Kasner exponents (A:i, A: 2 , ^ 3 ). The initial conditions are: p,iCi = 7 r/ 2 , 
Pi = 100, with i = 1, 2, 3 
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FIG. 8 : Evolution of expansion and volume. The initial conditions are: piCi = 7 r/ 2 , pi = 100, with 
i = 1,2,3. 


In figure we plot the density (p/pcrit) as a function of cosmic time. It shows how the 
density changes as the momentum of the held changes. One interesting thing to note is 
that the density at the bounce (t = 0 ) approaches the critical density (p/pcrit ~ 1 ) as the 
momentum of the held becomes larger, while the density of the bounce seems to increase as 
the momentum decreases. Recall that we are in the isotropic limit, so the change in density 
at the bounce can be interpreted as a transition from the k=l to the k=0 cases. Recall 
that the k=l FLRW model possesses an absolute bound on density only when inverse triad 
corrections are introduced, but the density at the bounce is always larger than in the k=0 
case where the density is always bounded by pcrit m- Additionally, the time between each 
bounce is longer as the momentum of the held increases, since the universe reaches larger 
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FIG. 9: Time evolution of the density (p/pcrit) for solutions with different values for the momentum 
of the field p^j,. The bounce is at t=0. The initial conditions for the three solutions are: fliCi = 7r/2, 
Pi = 1000, Pi = 10000, Pi = 100000, with i = 1,2, 3. Note that as the momentum p^ increases, the 
density at the bounce decreases and approaches the value of the k=0 FLRW model 


values of the volume and, therefore, it takes more time to recollapse. Thus, as one performs 
simulations with larger values of the momentum of the field one approaches, as a limiting 
case, the isotropic flat FLRW universe. This universe is not embedded within the Bianchi 
IX model but is rather a limit for the solutions, which can be approached as much as one 
wants in a hnite interval of time. The initial conditions for each of the solutions in Fig. 
are: /ijCj = n/2, pi = 1000, with ^ 4 x 10"^; /ijCj = 7i/2, pi = 10000, with p^ ^ 1 x 10® 
and fiiCi = 7 i/2 , pi = 100000, with p^ ^ 3 x 10^. 

Let us end this part by clarifying the situation at hand. Strictly speaking, we are not 
studying the isotropic limit corresponding to FLRW with k = 1, given that we choose 
exactly isotropic initial conditions, while in the k = 0 case, it is truly a limiting case since 
the FLRW model is not contained within the Bianchi IX model. Once we have cleared up 
what we mean by the isotropic limit, we proceed in the next section to study the solutions 
near to the isotropic limit, which we call the non-shear limit, since the shear has a small 
contribution throughout the dynamical evolution. 


D. Non-Shear Limit 

In the previous section we clarified that the isotropic limit is really made of isotropic solutions 
between the Bianchi IX, given that we chose exactly isotropic initial conditions. Now, we 
want to study the evolution of full Bianchi IX near to the isotropic solutions, namely, we 
want a small contribution from the shear. In order to do that, we choose initial data that 
contain exactly zero shear at the bounce, and then study how the shear changes during time 
evolution. While it is not a generic solution, it is an interesting limiting case. Furthermore, 
it can also be used to test the accuracy of the numerical code by checking whether the 
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FIG. 10: Evolution of volume and scale factors. The initial conditions are: /iici = /i 2 C 2 = MsCs = 
7r/2, Pi = 10000, p2 = 20000, ps = 30000. 


Shear Shear 




FIG. 11: During evolution the shear remains small and vanishes at the bounce due to the initial 
conditions: pici = / 22 C 2 = P 3 C 3 = 7 r/ 2 , pi = 10000, p 2 = 20000, p 3 = 30000. 


dynamics is symmetric across the bounce. The dynamics is expected to be time-symmetric 
given that all velocities {pi) vanish at the bounce point, and therefore there is no difference 
in time evolution forwards or backwards. We choose the initial conditions at the bounce for 
the effective equations and evolve them back and forward in time, the initial conditions are: 
PiCi = /i 2 C 2 = /I 3 C 3 = 7r/2, Pi = 10000, P 2 = 20000, ps = 30000. 

Figure [TO] shows the evolution of the volume and the scale factors. It is clear from the hgure 
that the evolution is symmetric across the bounce. This fact can be seen to validate the 
numerical code, but more intriguing, this solution is very different from the isotropic solution 
shown in the previous section. A new feature appears in the figure IT where it is shown the 
evolution of the shear, which vanishes at the bounce point, due to the initial conditions, and 
is different from zero during the evolution. The shear has a small value and the maxima are 
less than 4 x 10“^ (recall that the units are: c = 1, h = 1, G = 1,), which is a very small value. 
Note that the shear has four maxima and three minima, with one minimum being exactly 
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Volume 



FIG. 12: Time evolution of the volume for Bianchi IX and Bianchi I, effective and classical. Note 
that the solutions coincide at time t = 1835 where the initial conditions are chosen: /2iCi = 
6.0 X 10-^ Ji 2 C 2 = -1.1 X 10-^ liscs = 1.3 X 10-^ pi = 6.0 x 10^°, p 2 = 4.0 x 10^^, pa = 9.0 x lO^o. 


zero and there are two maxima with a larger value. If we consider initial conditions for 
PiiP 2 iP‘i with larger values, then the maxima of the shear become smaller than the current 
solution, shown in hgures 10 and 11 Next, we shall add a level of complexity and study the 
Bianchi I limit in what follows. 


E. Bianchi I Limit 

In order to hnd those Bianchi IX solutions that are close to the Bianchi I solutions we employ 
the following strategy: 

• Impose the initial conditions at the bounce for the effective Bianchi I equations, such 
that the held momentum is large. 

• Evolve forward in time the Bianchi I equations until the classical region. 

• Take the values from the Bianchi I evolution and use them as initial conditions for the 
Bianchi IX effective equations. 

• Evolve back in time the effective Bianchi IX equations. 

Using this procedure we get Bianchi IX solutions that start as a Bianchi I. The initial 
conditions that come from the Bianchi I evolution are: piCi = 6.0 x 10“®, Ji2C2 = —1.1 x 10“®, 
/i3C3 = 1.3 X 10“^, Pi = 6.0 X 10^°, p2 = 4.0 X 10^^, p3 = 9.0 x 10^°, the momentum of the 
held is p 0 ~ 1.83 x 10^^. The initial time is t = 1835 and the Bianchi IX equations are 
evolved back and forward in time. The momentum of the held for the Bianchi IX constraint 
is ~ 3.66 X 10^^, note that it is not exactly equal to the Bianchi I constraint, this is due 
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FIG. 13: Evolution of the Kasner exponents (ki) and the parameters that measure the dynamical 
contribution of matter (fi), anisotropies (S^) and intrinsic curvature (K). Recall that these quan¬ 
tities are not well defined at the bounce, so what we see is a transition from one value on one side 
of the bounce to a different value at the other side. The initial conditions at the time t = 1835 
are: flici = 6.0 x 10“^, //2C2 = —1.1 x 10“^, /23C3 = 1.3 x 10“^, pi = 6.0 x 10^^, p 2 = 4.0 x 10^^, 
P3 = 9.0 X 10^0. 
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to fact that Bianchi IX solutions should not necessarily follow the Bianchi I solutions the 
entire time, even in General Relativity. It is shown in hgure[^ where can be seen that the 
Bianchi IX and Bianchi I solutions coincide only at the initial conditions, and both follow 
the classical behaviour for each corresponding model. 

Figure [l^ depicts the evolution of the Kasner exponents and the parameters C, K, which 
remain almost constant in the region after the bounce. This is what we expect from a 
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FIG. 14: Time evolution of the curvature parameter [K) for the effective and classical Bianchi 
IX. Note that K is near to zero but has a non trivial behavior, as it becomes negative even at 
classical level. In the right figure we make a zoom up until the point where the effective and 
classical dynamics differ, near the bounce (recall that the quantity K is not well defined at the 
bounce). The initial conditions are: /iici = 6.0 x 10“®, li 2 C 2 = —1.1 x 10“^, /I 3 C 3 = 1.3 x 10“^, 
pi = 6.0 X 10^°, p2 = 4.0 X 10i\ p3 = 9.0 X 10^°. 
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FIG. 15: Evolution of the Xi quantities that contribute to the intrinsic curvature. The initial 
conditions at the time t = 1835 are: Ji\Ci = 6.0 x 10“^, JI2C2 = —1-1 x 10“®, /I 3 C 3 = 1.3 x 10“^, 
pi = 6.0 X 10^°, p2 = 4.0 X 10l^ p3 = 9.0 X 10^°. 


solution that evolves as a Bianchi I universe. The Kasner transitions in Bianchi I were 
studied in detail in |3H]- In figure 13 it can be seen than the dynamical contribution from 
the intrinsic curvature is near zero [K ~ 0), which is a condition that Bianchi I universes 
satisfy, but it has a non trivial behaviour, figure M Note that K becomes negative even at 
classical level. 

The Xi are the quantities that contribute to the intrinsic curvature Eq. (2.5), whose 
evolution is shown in figure I^, where it can be seen that ~ 0 after the bounce, while 
at the bounce Xi is large, and the other two are close to zero (x 2 ,X 3 ~ 0). Recall that 
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these are precisely the conditions that dehne the Bianchi II limit. Therefore, at the bounce 
there is a Bianchi II transition. The evolution before the bounce is highly non trivial, the 
dynamical contribution is changing constantly between T? and K (Fig. 13), now xi and X 2 
are contributing to K and X 3 remains close to zero, which we have classihed as the “Bianchi 
VIIq limit”. Of course, as we mentioned before, the Bianchi VIIq is just a label and needs to 
be properly studied to determine whether this dynamics really correspond to the effective 
dynamics of the loop quantum cosmology Bianchi VIIq model. 


To summarize, in Figs 
stages. At hrst the system evolves close to the Bianchi VIIq limit {x^ 


13 and it can be appreciated that the evolution has different 

0 and iF 7 ^ 0 ). 


Later on, it has a Bianchi II transition (x 2 , X 3 ~ 0 and xi 7 ^ 0) at the bounce and, hnally, it 
evolves like a the Bianchi I {K ~ 0, ki constants). This is a full Bianchi IX effective 

dynamics that can be characterized with some phases and transitions. Note that one can 
not assure that all the solutions can be characterized in this way; indeed it is highly probable 
that the dynamics in the full space of solutions would have even richer behaviours. 

Finally, we want to add some comments about the evolution of fl, T?, K in Fig. 13 First, the 
discontinuity of K is because they are ill-dehned at the bounce (they are dehned as a 

ratio over 9). Second, it is not completely clear in the hgure, but the curvature parameter K 
is negative in some interval of time before the bounce. Third, the evolution before t = —5000 
could not be performed due to the accumulation of numerical errors. Thus, it is not possible 
to try to reach any conclusions without a better numerical integration. This is due to the 
equations of motion being stiff and needing a special treatment for long time evolutions. 


V. CONCLUSIONS 


In this manuscript, we studied the numerical solutions to the effective Bianchi IX equations 
in loop quantum cosmology. This allowed us to answer some of the questions that motivated 
this work. Let us start by summarizing our results. Throughout this work, we considered 
two sets of effective equations, differing in their treatment of inverse triad corrections. We 
started with the question regarding the resolution of the big bang singularity. In sections 


IV A and IV B we showed that both effective theories in the large held momentum limit 


describe the same dynamics, resolve the singularity, and replace it by a bounce, and far 
from the bounce reproduce the classical dynamics. Furthermore, given that the Bianchi 
IX model has positive spatial curvature, then there is a inhnite number of bounces and 
recollapses. 


We found in sections IIV Cl and IIVEI that the set of effective Bianchi IX has as a limit the 
effective Bianchi I, and the effective isotropic FLRW with /c = 0. The FLRW with fc = 1 is 
contained within the Bianchi IX model only if the inverse triad corrections are not included 
(for the case with lapse N=V). When they are included (for lapse V=l), the FLRW is not 
embedded in Bianchi IX, but it is rather a limiting case when the volume at the bounce is 
large. 


Additionally, in section IV E we showed that the effective dynamics of Bianchi IX can be 
characterized by different phases. These phases can be classihed according to the values of 
the parameters x. 


introduced in section |IIB 
Moreover, the Kasner parameters k 


which contribute to the intrinsic curvature. 


are useful to characterize the Bianchi IX dynamics. 
In our case, the universe evolves from a universe close to Bianchi VIIq (X 3 ~ 0), later on 
it undergoes a bounce in which it exhibits a Bianchi II type transition (x 2 ,X 3 ~ 0), and 
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finally evolves close to a Bianclii I solution (xi, 0 : 2 , 2:3 ~ 0 and fci, A; 2 , ^3 constants) far from 
the bounce. Note that this is the hrst time that the Bianchi VIIq model is mentioned in 
the context of loop quantum cosmology. This model has not been studied yet, but we 
are speculating that the solutions that we get are close to what would be the effective 
dynamics of Bianchi VIIq, which would come from the “improved dynamics” Bianchi VIIq 
loop quantization. 

All the result that we have obtained are in the limit of large held momentum, where both 
effective theories of Bianchi IX have the same dynamics. If these two theories really describe 
the quantum dynamics of the semiclassical states, as it is expected, then our results indicate 
that the semiclassical dynamics of both quantum theories describe the same physics for a 
realistic large universe. 

We shall end by discussing some related issues to what we have studied in this manuscript. 

1. When one considers many bounces and recollapses, then new features and interesting 
behaviours appear, when looking at the dynamics of the observables [28], especially in 
the vacuum case that is probably the most interesting one. 

2 . One can also perform qualitative studies that are independent of the matter content 
of the theory, by just imposing energy conditions. In that case, one can explore 
what happens to the BKL behaviour when quantum corrections -as captured by the 
effective dynamics- are incorporated. Those rather intriguing results will be published 
in a forthcoming communication [27] . 

3. In this manuscript we have considered only two effective theories of Bianchi IX, but 
there are more effective theories [IllET] which explore the ambiguities in the quanti¬ 
zation of Bianchi IX. 

4. Preliminary studies 1271 indicate that our results will also apply to the other effective 
theories [THEI] in the limit of large held momentum. 
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Appendix A: Effective Equations of Motion 

In this section we collect all the effective equations of motion that we use in the numerical 
study of the solutions. Some of these equations of motion come directly from the original 
articles, we just want to put all together in this appendix. The programs developed are 
available upon request. 
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Flat FLRW, k = 0 


The effective Hamiltonian is [21 ES] 


Ch 


k=0 


3 


sin2(A/?) 



(Al) 


where the lapse is iV = H, with V the physical volume of the fiducial cell. The Poisson 
brackets are {/S, V} = AtiGj and {0,^,/,} = 1. The variables (V, (3) are related with (c, p) by 


1 / = ^ p 


The equations of motion for (V,/3) are 


c 

Vp ’ 


V = sm{X/3) cos(A/3), 

yA 

P = —^ H sm\XP). 


(A2) 


(A3) 

(A4) 


Closed FLRW, k = 1 


For the closed FLRW model there are different quantum theories that depend on the quan¬ 
tization method used for defining the curvature [7], and also depending on whether the 
inverse triads corrections are included or not. For our study we choose the equations that 
come from the quantization where the curvature is calculated using the connection, which 
is obtained from the ‘open’ holonomies. Even so the freedom to choose the inverse triad 
correction remains. Therefore we are going to consider two effective theories, depending of 
the lapse function. The effective Hamiltonian with lapse N = V is given by [a ED] 


C 


( 1 ) _ 


Hk 


[sin" A^) - 2 Dsm\D + (1 + 7")^"] -J ~0 


The equations of motion are 


V = —V cos{X/3) [sin(A/?) — D ], 
yA 
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^ ^ \ „ [3 sin^ A/3 — AD sin A/j -|- (1 -|- y^)T>^1 — 27rG’y ^ . 

o.. \ ^ L J yz 


2yA2 

The effective Hamiltonian with lapse iV = 1 is HlEol 


C 


S'- ^ A/3 + (1 + - y ~ 0 , 


with 


A(P) = yy(R + W-|R-W|) = 


r v/Vc , iiv<v, 

11 ,ifR>Rc 


(AS) 

(A6) 

(A7) 

(AS) 

(A9) 
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and Vc = 27r7A£pi. The equations of motion are 
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cos A/?(i/sin A/? - , 

7 A 


(AlO) 


27 A 2 


VAy + A(i/)^ sin^ A/3 - 2A'd sin A/3 (^A{y)V-^/^ + AyV^/^^ 


+AH\l + f 
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A,y = :rT^ ( 1 


2i4 


|V^-K 

i/-i4 


/ 1 /K , ifi/< K 

\0 ,ifV^>K 


(All) 


(A 12 ) 


In both theories the dehnition of the (V, /3) variables is the same than the flat FLRW case 
iy = = c/^yp, {/3,V} = 47 rG 7 ), but now the volume correspond to the physical 

volume of all the space. 


Bianchi I and II 

The effective Hamiltonians for Bianchi I and II can be written in one single expression 
P 1 P 2 P 3 


^Hbii 


87rG72A2 L 
1 


sm piCi sm P 2 C 2 + sm P 2 C 2 sm + sm sm piCi 


a(P2P3)^A 


n , 2^/«P2P3V 


_^„0, (A13) 


where a allows to choose between Bianchi I (a = 0) and Bianchi II {a = 1). The Poisson 
brackets are = SiiG'ySj and {</>,P 0 } = 1. The equations of motion are |32] 


Pi = 

T)^ 

—— (sin P2C2 + sin P3C3 + Ax) cos piCi, 

(A14) 
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P2 = 

^ (sin/iiCi + sin/iaCa) COSP2C2, 

7/^2 

(A15) 

P 3 = 

Jp‘ 

^ (sin piCi + sin P2C2) cos P3C3, 

(A16) 
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P 2 P 3 
' 27 A 2 L 


2(sin piCi sin P2C2 + sin piCi sin P3C3 + sin P2C2 sin P3C3) 


+P 1 C 1 cos piCi(sin P 2 C 2 + sin P 3 C 3 ) - P 2 C 2 cos p 2 C 2 (sin piCi +sinp 3 C 3 ) 

-P 3 C 3 cosp 3 C 3 {smpiCi + smp 2 C 2 ) + A^a;^(l + 7 ^) 

+Ax{piCiCos piCi — sin piCi) , (A17) 
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C 2 = 


PlPS 

'27A2 


2 (sin fiiCi sin/i 2 C 2 + sin fXiCi sin + sin/i 2 C 2 sin P 3 C 3 ) 


-piCi cos/iiCi(sin/i2C2 + sin/iaCa) + P2C2 cos/i2C2(sin/iiCi +sin/i 3 Ca) 

-P3C3 cos/2aCa(sin/iiCi + sin/i2C2)] - + 7^) 

—Ax(/iiCi cos/iiCi — 3 sin/ 2 iCi) , (A 18 ) 


C3 = 


with 


P1P2 

'27A2 


2 (sin /iiCi sin /i 2 C 2 + sin /xiCi sin /xaCa + sin fi 2 C 2 sin /xaCa) 


-/iiCi cos/iiCi(sin/i2C2 + sin/X3Ca) - P2C2 cos/i2C2(sin/iiCi + sin/i3C3) 
+/i3C3 cos/i3C3(sin/iiCi + sin/i2C2)] - A^a;^(l + 7^) 

—Aa;(/iiCi cos/iiCi — 3 sin piCi) 


X = a 


P2P3 

Pi 


(A19) 


(A20) 


Bianchi IX 


The effective Hamiltonian for Bianchi IX with lapse A^ = H is [18 
P1P2P3 


^Hbix 


87 rG 7 ^A^ 

({P\P2)^^‘^ 


( sin /iiCi sin /i 2 C 2 + sin /i 2 C 2 sin paCa + sin P 3 C 3 sin piCi) 


47 rG' 72 A V 


. - , iP2P3f^‘^ . - , {P3Plf^'^ . - 

sin /xaCa H-—— sin /xiCi H-—— sin /r2C2 


i^^(l + 7 ") 
87 rG 7 ^ 

4«. 


The eqnations of motion are [18 
1 


‘^{pI+pI + pD - 


y/pl 

P1P2 

P3 


VP2 


P2P3 

Pi 


P3P1 

P2 


2 1 


Pi = 


Cl = 


7 


^(sinp 2 C 2 + sinpaca) + 2 i 9 p 2 P 3 

/ii 


cos /iiCi 


(A21) 


(A22) 


1 /P 2 P 3 , 

■ - < —sm /ilCi sm /i2C2 + sm piCi sm P3C3 + sm P2C2 sm P3C3 

7 I V 

, hlCl _ , . _ . _ X /I2C2 _ , . _ _ X 

cos piCifsm /i2C2 + sm P3C3) -— cos /i2C2(,sm /xiCi + sm /iaCa) 


2 

P3C3 


2 cos/i 3 C 3 (sin/xiCi + sin/i 2 C 2 )j +2191^ — 


P1P2 . _ , P1P3 . - 

sm/xaCa H-sm/i 2 C 2 


P3 


P2 


P2P3 

3pi 


sm fiiCi 


lP 2 P 3 


2 Pi 


1 

2^ 


+ -^-^cicospici - -P2C3 cos paCa - -P3C2 cos/i2C2 I 


1 

2^ 


+i 9'(1 + 7 ') 


P2 , pi 


4 pi — 2 pi “2 H 2 J ^ 3 

pi P 2 J Pi J 


,pIpI 


(A23) 
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The other equations of motion are obtain by permutations. The effective Hamiltonian for 
Bianchi IX with lapse iV = 1 is given by [HI [30] 


C 


( 2 ) 


V^A{V)h\V) 

87rG72A2I/6 

dA{y)h\V) 


(sin /iiCi sin /i 2 C 2 + sin piCi sin psCs + sin P 2 C 2 sin P 3 C 3 ) 


.,22, 


,2„2 , 


AttG^^XV^ (^PiP 2 sin + pIpI sin piCi + plpl sin p2C^ 

d^{l + -f^)A{V)h'^{V) 


SvrGy^I/'i 


X (^v{ pI +pI+ pI) - [{pmT + (piPsT + iP2P3V] 


h%v)V^ 

Wf~ 


pI-o, 


with 


h{v) = ^yv+v, - 

The derivatives of h{V) and AiV) are given by 


j_ i PjPk 

0 , if ^ > K 


= { ^-\ Pi . ,iyjykyi 


and 


1 IpjPk 

= 4 V ^ 


1 

^/vTVc V -V, 


3 ykyi 


The equations of motion are HHIEQ] 


1 _ \v^h^{y)pi 

p'l = —cospiCi - (sin/i2C2 + smpsCs) - 2'dp2Pz 






(A24) 

(A25) 

(A26) 

(A27) 


(A28) 
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Cl = - 


h^{V) 

7A2\/6 


2piplpiA{y)h{y) + V^A^p^h{V) + QV^A{V)h 


.pi 


X (sin piCi sin P2C2 + sin piCi sin + sin P2C2 sin P3C3) 


2i9 


h%V) 


{2piplA{V)h{V) + plplA^p^h{V) + 4plplA{V)h^p^) sin/igCs 


{2piplA{V)h{V) + plplA^p^h{V) + AplplA{V)h^p^) sin/isCa 
{P2PI^,pMV) + 4P2P3^(^)^,Pi) sin piCi 


3/2 3 / 2 ' 
P 2 Ps 


- y V'^^P (sin/i2C2 + sin/igCg) - 




C2COSP2C2 f pl^'^pl^'^pl^'^h'^iV) 

2v;^ I 


(sin/iici + sin/igcs) - 2§p3V^ 


3lTn;,4^Tn^3COS/i3C3 . „ , • - 1 OQ T/ 

Al(i/)/i [V)^— - ---(sm/iiCi + sin/iaCa) - 2i9p2V" 


2V;'^7 




i9^(l + 7 ^) 


h%V) 


4p,A{V)Vh{V) - —pih*^{V)A{V){pt + pt) 


(pI + pI+pD 


P 2 P 3 

Pi 


A{V)h{V) + 8A{V)Vh,p^ + 2A^p^Vh{V) 


- ^ ph\V)h^p^A{V) + h\V)A^p^ (^p^^p^ + p^pl + p^pl^ 

The other equations of motion are obtain by permutations. 


t^GipI 

PWP 1 P 2 P 3 


(A29) 


Appendix B: Convergence and Conservation 
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FIG. 16: Convergence of the Hamiltonian constraint and field momentum conservation. The initial 
conditions are: pici = 7r/6, P 2 C 2 = 2tt/5, P 3 C 3 = Svr/S, pi = 12000, p 2 = 16000, p 3 = 14000. 













FIG. 17: The first plot shows the error functions Li and L 2 . The second plot shows the convergence 
order for three different samples. The initial conditions are: fiici = 7 r/ 6 , fi 2 C 2 = 27r/5, /U 3 C 3 = hvr/S, 
Pi = 12000, p2 = 16000, p 3 = 14000. 


There are two important points that need to be addressed in a numerical work: one is the 
convergence of solutions and the other is the evolution of conserved quantities. Numerical 
solutions must also evolve on the constraint surface and they must preserve the conserved 
quantities. This ensures that they are evolving on the physical phase space. In hgure IT 
we show the convergence of the Hamiltonian constraint {Ch ~ 0) and the held momentum 
evolution, which shows that is conserved. The quantities plotted in Figure 17 are the 
error functions 


Li = Meix\CH{t)2-CH{t)i\, L 2 = Max^|C|f(t)2 -C]j{t)i \, (Bl) 

where the sub-indexes in the Hamiltonian constraint mean resolution 2 (with dt 2 ) and res¬ 
olution 1 (with dfi), such that dt 2 = dti/2. The method used to integrate the equa¬ 
tions is a Runge-Kutta 4 (RK4), while the resolutions used for the convergence tests are 
dt = 0.032, 0.016, 0.008, 0.004, 0.002, 0.001. We can dehne the convergence order as 

n(.t) = (B2) 

with fi any evolved function at resolution i, with dti > dti+i. The convergence factor n for 
a RK4 must be u = 2^^ = 16 and we get in our solutions (n) ~ 16.6, where {n) denotes the 
average in time and taking the mean of different samples. 
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